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Abstract 
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monopole and core-induced contributions. The asymptotic behavior of the core-induced 
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values of the solid angle corresponding to strong gravitational fields. As an application of 
general results the flower-pot model for the monopole's core is considered and the expectation 
values inside the core are evaluated. 
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1 Introduction 



Symmetry braking phase transitions in the early universe have several cosmological consequences 
and provide an important link between particle physics and cosmology. In particular, different 
types of topological objects may have been formed by the vacuum phase transitions after Planck 
time [H[2]- These include domain walls, cosmic strings and monopoles. A global monopole is 
a spherical symmetric gravitational topological defect created by a phase transition of a system 
comprised by self-coupling scalar field, ip a , whose original global 0(3) symmetry is sponta- 
neously broken to U(l). The matter fields play the role of an order parameter which outside the 
monopole's core acquires a non-vanishing value. The global monopole was first introduced by 
Sokolov and Starobinsky [3J. A few years later, the gravitational effects associated with a global 
monopole have been considered in Ref. [4], where the authors have found that for points far 
from the monopole's center, the geometry is similar to the black-hole with a solid angle deficit. 
Neglecting the mass term we get the point-like global monopole spacetime with the metric tensor 
given by the following line element 

ds 2 = dt 2 - dr 2 - a 2 r 2 (d6 2 + sin 2 6d(j) 2 ) , (1) 

where the parameter a 2 is smaller than unity and depends on the energy scale where the sym- 
metry is broken. It is of interest to note that the effective metric produced in superfluid 3 He — A 
by a monopole is described by the line element (pQ) with the negative angle deficit, a > 1, which 
corresponds to the negative mass of the topological object [5]. The quasiparticles in this model 
are chiral and massless fermions. 

In quantum field theory the non-trivial topology of the global monopole spacetime induces 
non-zero vacuum expectation values for physical observables. The quantum effects due to the 
point-like global monopole spacetime on the matter fields have been considered for massless 
scalar [6j and fermionic [7J fields, respectively. The influence of the non-zero temperature on these 
polarization effects has been discussed in Ref. [8] for scalar and fermionic fields. Moreover, the 
calculation of quantum effects on massless scalar field in a higher dimensional global monopole 
spacetime has also been developed in Ref. [9]. The quantum effects of a scalar field induced 
by a composite topological defect consisting a cosmic string on a p-dimensional brane and a 
(m + l)-dimensional global monopole in the transverse extra dimensions are investigated in Ref. 
|10j . The combined vacuum polarization effects by the non-trivial geometry of a global monopole 
and boundary conditions imposed on the matter fields are investigated as well. In this direction, 
the total Casimir energy associated with massive scalar field inside a spherical region in the 
global monopole background has been analyzed in Refs. |11[ [12] by using the zeta function 
regularization procedure. Scalar Casimir densities induced by spherical boundaries have been 
calculated in Refs. |X3|, [T4] to higher dimensional global monopole spacetime by making use of 
the generalized Abel-Plana summation formula [15]. More recently, using also this formalism, 
a similar analysis for spinor fields obeying MIT bag boundary conditions has been developed in 
Refs. PEIII3. 

In general, the quantum analysis of matter fields in global monopole spacetime consider this 
object as been a point-like one. Because of this fact the renormalized vacuum expectation value 
of the energy-momentum tensor presents singularities at the monopole's center. Of course, 
this kind of problem cannot be expected in a realistic model. So a procedure to cure this 
divergence is to consider the global monopole as having a non-trivial inner structure. In fact, 
in a realistic point of view, the global monopoles have a characteristic core radius determined 
by the symmetry braking energy scale where the symmetry of the system is spontaneously 
broken. A simplified model for the monopole core is presented in Ref. [18]. In this model 
the region inside the core is described by the de Sitter geometry. The vacuum polarization 
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effects associated with a massless scalar field in the region outside the core of this model are 
investigated in Ref. [19]. In particular, it has been shown that long-range effects can take place 
due to the non-trivial core structure. Recently the quantum analysis of a scalar field in a higher- 
dimensional global monopole spacetime considering a general spherically symmetric model for 
the core, has been considered in Ref. [20] . In the four-dimensional version of this model the 
induced electrostatic self-energy and self- force for a charged particle are investigated in Ref. [21] . 
Continuing in this direction, in the present paper we analyze the effects of global monopole core 
on properties of the fermionic quantum vacuum. The most important quantities characterizing 
these properties are the vacuum expectation values of the energy-momentum tensor and the 
fermionic condensate. Though the corresponding operators are local, due to the global nature of 
the vacuum, the vacuum expectation values describe the global properties of the bulk and carry 
an important information about the structure of the defect core. In addition to describing the 
physical structure of the quantum field at a given point, the energy- momentum tensor acts as 
the source of gravity in the Einstein equations. It therefore plays an important role in modelling 
a self-consistent dynamics involving the gravitational field. The problem under consideration is 
also of separate interest as an example with gravitational polarization of the fermionic vacuum, 
where all calculations can be performed in a closed form. The corresponding results specify the 
conditions under which we can ignore the details of the interior structure and approximate the 
effect of the global monopole by the idealized model. The exactly solvable models of this type 
are useful for testing the validity of various approximations used to deal with more complicated 
geometries, in particular in black-hole spacetimes. 

The plan of this paper is as follows. Section [2] presents the geometry of the problem and the 
eigenfunctions for a massive spinor field. By using these eigenfunctions, in Section[3]we evaluate 
the vacuum expectation values of the energy-momentum tensor and the fermionic condensate 
in the region outside the monopole core. In Section 0] we consider the special case of the flower- 
pot model for the core and the vacuum expectation values are investigated in both exterior 
and interior regions. Various limiting cases are considered. The main results of the paper 
are summarized in Section [5j In Appendix [A] we discuss the contribution of possible bound 
states into the vacuum expectation values of the energy-momentum tensor and show that in 
the flower-pot model there are no bound states. Throughout we use the system of units with 
h = c=l. 

2 Model and the eigenfunctions for the spinor field 

In this section we analyze the eigenfunctions for a massive fermionic field in background of the 
global monopole geometry considering a nontrivial inner structure to the latter. The explicit 
expression for the metric tensor in the region inside the monopole core is unknown and here 
we consider a general model for a four-dimensional global monopole with a core of radius a, 
assuming that the geometry of the spacetime is described by two distinct metric tensors in the 
regions inside and outside the core. Adopting this model we can learn under which conditions 
we can ignore the specific details for its core and consider the monopole as being a point-like 
object. In spherical coordinates we will consider the corresponding line element in the interior 
region, r < a, with the form 

ds 2 = e 2u(r Ut 2 - e 2v ^dr 2 - e 2h ^\d9 2 + sin 2 fl# 2 ) . (2) 

In the region outside, r > a, the metric tensor is given by the line element ([T|), where the 
parameter a codifies the presence of the global monopole. For an idealized point-like global 
monopole the geometry is described by line element ([1]) for all values of the radial coordinate 
and it has singularity at the origin. 
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In Eq. ([5]) the functions u(r), v(r), h(r) are continuous at the core boundary, consequently 
they satisfy the conditions 

u{a) = v(a) = 0, h(a) = ln(cra) . (3) 

If there is no surface energy-momentum tensor on the bounding surface r = a, the radial 
derivatives of these functions are continuous as well. When the surface energy-momentum 
tensor is present the junctions in the radial derivatives of the components of the metric tensor 
are expressed in terms of the surface energy density and stresseqj. Introducing a new radial 
coordinate f = e h ^ with the core center at r = 0, the angular part of line element ([2]) is 
written in the standard Minkowskian form. However, with this choice, in general, we will 
obtain non-standard form of the angular part in the exterior line element ([I]). In the model 
under consideration we will assume that inside the core the spacetime geometry is regular. In 
particular, from the regularity of the interior geometry at the core center one has the conditions 
u(r),v(r) — > and h(r) ~ lnf for f — * 0. 

We are interested in quantum effects of a spinor field propagating on background of the 
spacetime described by line elements ([TJ and ([2]). The dynamics of the massive spinor field in 
curved spacetime is governed by the Dirac equation 

i'fV^ -Mip = , (4) 

with the covariant derivative operator 

v M = + r M . (5) 

Here 7^ = e^7^ are the Dirac matrices in curved spacetime, and is the spin connection 
given in terms of the flat 7 matrices by the relation 

= \l {a) l {b) e\ a) e {b) ^ . (6) 

In the equations above ej^ is the vierbein satisfying the condition el L\ e ^V = g^ v . 
In the region inside the monopole core we choose the basis tetrad 



u 

(a) 



/ e~ u ^ \ 

e ~ v ( r ) gin 9 cos ip e~ h ^ cos 8 cos <p —e~ h ^sm(p/sm6 

e~ v ( r ^ sin 9 sin ip e~ h ^ cos 9 sin ip e~ h ^ cos <p/ sin 9 

\ e~ v ^cos8 -e- h(r hin9 / 



(7) 



where the rows of the matrix are specified by the index a and the columns by the index \i. By 
using Eq. (|7|), for the non- vanishing components of the spin connection we find 

To = \u'e u - v 7 (0) 7-r, 

T 2 = l -(l-h'e h - v )t-0, (8) 
T 3 = - l -sm9{l-h'e h - v )t-9 , 

where the prime denotes derivative with respect to the radial coordinate, f , 9 and (p are the 
standard unit vectors along the three spatial directions in spherical coordinates, and 



a 
a 



(9) 



1 This point will be analyzed in more detail at the end of the next section 
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with a = (o"i, o"2, 03) being the Pauli matrices. From the obtained results, for the combination 
entering in the Dirac equation we have 

rr^ = ^^7-r-e- h (l-h'e h - v )^-f . (10) 
After some intermediate steps, the Dirac equation reads 

ie-V 0) d t V + ie- u/2 - v - h 7 (r) d r (e u/2+h ip) - ie~S {r) (S ■ L + l)tp — Mip = , (11) 

with L being the standard angular momentum operator. 

To find the vacuum expectation value (VEV) of the energy-momentum tensor we need the 
corresponding eigenfunctions. For the problem under consideration there are two types of 
eigenfunctions with different parities which we will distinguish by the index a = 0,1. These 
functions are specified by the total angular momentum j = 1/2, 3/2, ... and its projection 
m = —j, — j + 1, j. Assuming the time dependence in the form e~' lujt , the four-component 
spinor fields specified by the set of quantum numbers (3 = (akjm) with k 2 = uj 2 — M 2 , can be 
written in terms of two-component ones as 

where 

n a = (-iy, i« = j~Y> & = i + T' (13) 

and £lji aUl (6,(p) are the standard spinor spherical harmonics (see, for instance, Ref. |22j). The 
latter are eigenfunctions of the operator fC = a ■ L + / as shown below: 

KSljlam = -KvQjlrm , K a = -Tl a (j + 1/2). (14) 

Note that we have the relation 

^ji' a m = i la ~ l '°{n ■ a)nji„ m ■ (15) 

Substituting the function ipp into the Dirac equation above, and using for the flat Dirac 
matrices the representation given in Ref. [22], for the radial functions we obtain a set of two 
coupled linear differential equations: 

f p + (u>/ 2 + h' + K a e v ~ h )fp = {M + e- u ui)e v gp , (16) 
^ + ( u '/2 + h'- Ka e v - h )gp = (M-e- u iv)e v fp . (17) 



In the region r > a for the radial functions we have the solutions 

Mr) = fj 3 Cx) (r) = ^=[c 1 J, a (kr)+c 2 Y^(kr)], 



(18) 



where J Va {x) and Y Va (x) are the Bessel and Neumann functions of the order 

j + 1/2 n a 

v a = 7C- 

a 2 



(20) 
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By taking into account relation (|15|) . the corresponding eigenfunctions are written in the form 

-,(«) _ e -iut ( f^( r )^jl,m(0, <p) 



^ x - e \ . m / :/: j z"? ' r L , ) • (2i) 



In Eqs. (|18p and (|19p . the integration constants Cx and C2 are determined from the matching 
condition with the interior solution. 

The regular solution to Eqs. (|16p . (117j) in the interior region, r < a, we will denote by 



//9(r) = iWnfc), 9p(r) = --^R 2 , ntr ( r ,k). (22) 

Near the core center, f — ► 0, these functions behave like i?i,n CT ^ r y+ ( 1_n<T -' //2 and -R2,n CT 
fi+( 1 +™<r)/2^ wriere ^j^g rac li a i coordinate f is introduced in the paragraph after formula ([3]). 
Note that the functions Ri :TlcT (r,k) and n a kR2- na (r,—k)/(M — to) are solutions of the same 
equation and, hence, 

R%-na ( r > ~ k ) = const • Rl,n* ( r > k ). ( 23 ) 
The interior eigenfunctions have the form 

./.(in) _ -fart f Rl,nAr,k)n jltTm (0,ip) \ 

From the continuity condition of the eigenfunctions on the surface r = a, for the coefficients c\ 
and C2 in Eq. (|18|) one finds 

ci = -^n CT A;a 3/2 [i?i inCT (a,fc)Y" I/(T+n(T (A;a) -i? 2 ,n CT (a,A;)l^ (T (A;a)] , (25) 

C2 = ^n a ka 3/2 [R ljna (a,k)J l/a+ntT (ka) - R 2 ,n a (a,k)J u<T (ka)] . (26) 



Note that from Eqs. (|16() . (|17p . for the values of the interior radial functions on the boundary 
of the core we have the following relations 

- n CT fci? 2 ,n CT (a, k) = R[ n<T (a, k) + (u' a /2 + h' a - n a X/a) Ri, n „(a, k), (27) 
n (T kR 1<niT (a,k) = R 2n<j (a,k) + (u' a /2 + h' a + n a X/a) R 2 ,n a (a, k), (28) 

where Rj n ^(a,k) = dRj >n<7 (r,k)/dr\ r=a , u' a = du/dr\ r=a , h' a = dh/dr\ r=a , and we have intro- 
duced the notation 

A = (i + l/2)/a. (29) 

Substituting the expressions for the coefficients c\ and c 2 into the formulae for the radial 
eigenfunctions in the exterior region, one finds 

/ (cx) 

_ irR 1>n<J {a, k)e lU)t ( 

5^(7, o(^^) kv)£lji m \ fqr\\ 

^ 2^Ta V ^M9u„nAka, kr){n ■ a)n jlm J ' W 

with the notation 

9vv, P (x,y) = Y Uiy (x)J V{J+p {y) - JuA x ) Y vv+ P {y)i P = O,^. (31) 

Here and in what follows, for a function F V(T (x) with F = J, Y we use the tilded notation defined 
by the formula 

F Ua (ka) = -n a ka [F Va+na (ka) - F Ua (ka)R 2 ,n a (a, k)/R 1)Ua (a, k)] 

R 'l,nM k ) , 1 / , ./ 1 



kaF'ika) 



fl P 1 h\ + OU « + aH a ~ O 

•Rl.n^COjfc) 2 2 



i^fca). (32) 
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In deriving the second form we have used the relation (|27|) and the recurrence formula for the 
cylindrical functions. 

The eigenfunctions are normalized by the condition 



(33) 



where 7 is the determinant for the spatial metric and Spp is understood as the Kronecker delta 
symbol for the discrete components of the collective index (3 and as the Dirac delta function for 
the continuous ones. As the normalization integral is divergent for (5 = /3' , the main contribution 
comes from large values r. We may replace the cylindrical functions with the argument kr in 
Eq. (|30p by the corresponding asymptotic expressions. In this way one finds 



Rl,nA a ' k ) 



2k(M + uj) 
7r 2 a 2 auj 



Ji(ka)+Yl(ka) 



(34) 



This relation determines the normalization coefficient for the interior region. In addition to the 
modes with real k, modes with purely imaginary k can be present. These modes correspond to 
the bound states. The eigenfunctions for the bound states and their normalization are discussed 
in Appendix 1X1 



3 Vacuum expectation values in the exterior region 

In this section we consider the VEVs for the energy-momentum tensor and the fermionic con- 
densate in the region outside the global monopole core. We expand the field operator in terms 
of the complete set of eigenfunctions {ijA 1 , ^} : 

(3 

where fig is the annihilation operator for particles, and is the creation operator for antiparti- 

cles, ijjp = ipp, for uj > and ipp ^ = ipp, for u> < 0. In order to find the VEV for the operator 
of the energy-momentum tensor, we substitute the expansion (|35p and the analog expansion for 
the operator tp into the corresponding expression for spinor fields, 

t^^} = l -[h{^ v) i> - {V(»i>hu)i>] ■ (36) 

By making use of the standard anticommutation relations for the annihilation and creation 
operators, for the VEV one finds the following mode-sum formula 

(0|V|0) =J2T^{^\x),^-\x)} , (37) 

13 

where |0) is the amplitude for the corresponding vacuum state. 

Substituting the eigenfunctions ([30]) into the mode-sum formula ([37]) . using the formula 

£ |%z CTm (^)| 2 = (38) 

m=-j 
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and relation (|34|) . the VEV of the energy-momentum tensor is presented in the form (no sum- 
mation over /j,) 



mo) 



87ra 2 a 3 i 



i=i/2 

In formula (1391) we use the notations 



,_1 In (7=0,1 J0 



/ill [aJ,^j»(«.y)] 



Jl(x) + Yl(x) 
(y / x 2 + M 2 a 2 - Ma)g*{x, y) 



(39) 



+(Vx 2 + M 2 a 2 + Ma)^„ (x, y) 



(40) 



x 3 (2i/ (T + n CT ) 



[i&.ofa y) + Ot^fa V)] ~ 2 /<Sr [X, 9v a , V {x, V)\ , (41) 

9u a fi{x,y)gu a ,n a {x,y), fi = 2,3, (42) 



2yVx 2 + M 2 a 2 ' 

and the function gv a , P (x,y) is defined by Eq. (f3"Tj) . The expression on the right of Eq. (f39l) is 
divergent. It may be regularized introducing a cutoff function ip^u) with the cutting parameter 
\i which makes the divergent expressions finite and satisfies the condition — ¥ 1 f° r A 4 — ^ 0- 
After the renormalization the cutoff function is removed by taking the limit fi —* 0. In the 
discussion below we will implicitly assume that the corresponding expressions are regularized. 
The parts in the VEVs induced by the non-trivial structure of the core are finite and do not 
depend on the regularization scheme used. 

To find the part in the VEV of the energy-momentum tensor induced by the non-trivial core 
structure, we subtract the corresponding components for the point-like monopole geometry. The 
latter are given by the expressions which are obtained from Eq. (139() replacing the integrand by 
faul [x, J UtT (xr /a)] (see Ref. PS]). In order to evaluate the corresponding difference we use the 
relation 

fatl [x,JuAxr/a)\ = -- 2^ ~~j~sS~, J fvl x,HlJ{xr/a) 



J 2 (x)+Y 2 (x) 



Z2Hil\x) 



(43) 



where Hp (z), s = 1,2, are the Hankel functions. This allows to present the vacuum energy- 
momentum tensor components in the form 



(o|t:|o) = (o m |r;|o m ) + <t. 



)JL /C ) 



(44) 



where (0 m |T^|0 m ) is the VEV for the point-like monopole and the part (no summation over [/,) 



167TQ 2 a 3 r 



E^ + uE E 



3=1/2 



CT =0,1 8=1,2 



dx 



Hll\x 



■f^\x,H u lXxr/a) 



(45) 



is induced by the non-trivial core structure. In formula (|45p . the integrand of the s = 1 (s = 2) 
term exponentially decreases in the upper (lower) half of the complex plane x. Consequently, 
we rotate the integration contour on the right of this formula by the angle tt /2 for s = 1 and by 
the angle —ir/2 for s = 2. This leads to the representation 



(T't 



fj,)c 



167ra 2 a 3 r 



£(2j + i) Y, I dx^n 



i=i/2 
J u „(e^ 2 x) 



<T=0,1 



8=1,2 



/« 



3^ ni / 2 x,H u s \e r ' sWi/2 xr/c 



(46) 
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where rj s = (— l) s+1 . 

First of all let us consider the part of the integral over the interval (0, Ma). By using the 
standard properties of the Hankel functions it can be easily seen that in this interval one has 



inX f M \ e ^x,HjS^ 2 xr/a)\ = e~™ x fM \e~^ 2 x, (e^xr/a) 



(47) 



Further, from equations (|16p and (|17p it follows that the interior solution can be written as 
Ri,n a { r , k) = const • R(r,tu). On the base of this observation for the combination entering into 
the tilted notations in Eq. (|46p one finds 

R' hn ^a,e m / 2 x/a) R' hria (a, e'^x/a) _ _ 



Ri,n a (a, e™/ 2 x/a) R ljUa (a, e~™/ 2 x/a) 
Now, it can be seen that for these values x one has 

—in\ J»Ae m/2 x) j7rA Jvje-™ l 2 x) 



, < x < Ma. 



(49) 



Combining relations (I47h and (149p . we see that the part of the integral in Eq. (|46p over the 
interval (0, Ma) vanishes. 

To simplify the part of the integral over the interval (Ma, oo), we note that for the functions 
with different parities the following relation takes place: 

-e-^ A /2 , [e- n/2 x,^)(e- 1! / 2 ir/a)| , (50) 



m \e m l 2 x,H$(e^ 2 xr/a) 



where a, a' = 0, 1, and a ^ a'. Further we note that the functions R2,n tT ,( r , —k) and Ri t n a ( r , k) 
satisfy the same equation and, hence, i?2,n / (r, —k) = const • i?i, n£r (r, k). By using relations (|27|) 
and (|28|) . now it can be seen that 



R 2 ^,(a,ke^ 2 ) _ fl 1|W >,fce±™/ 2 ) 
i? li7V (a, ke^/ 2 ) R 2 , na (a, ke±™l 2 ) ' 
From this formula, by taking into account that zv = v a + n u, we find the relation 



fl-W( e «/2x) 



H {2) ,{e-™/ 2 x) 



(51) 



(52) 



Combining formulae (|50p and (j52|) . we see that the different parities give the same contribution 
to the VEV of the energy- momentum tensor. 

By taking into account this result and introducing the modified Bessel functions, for the 
core-induced part in the VEV we find (no summation over //) 



su oo 



2ir 2 a 2 i 



dx 



M 



l l/a-l/2\ ax ) ^fan.) 



Here for a given function f v {y) we use the notations 

' M 2 



s=X,2 K l/a~l/2( aX ) 

ir] s M 



[x,K, 



l/a-l/2 



xr 



(53) 



1 + 



1 



irj s M 



Vx 2 - M 2 



Vx 2 - M 2 

/, 2 +i(y) 



/ 2 (y) 



fiiv)-fi+i(v)-*f 



2u + l 

y 



fu(y)fu+i(y), 



x f 2u J r 1 fu(y)fu+i(y), m = 2 ,3, 

2y 



(54) 
(55) 
(56) 
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with Xk = —1 (the function F^' 713 ' [x,I u (y)\ with Xj = 1 will be used below), and the barred 
notation is defined by the formula 



f(^\x) =xf{x) 



R' (a,e^ m / 2 x/a) 

l l ' na ) rp: — -r-r + du' 2 + ah' a - 1/2 



f{x). (57) 



It can be checked that the core-induced part in the VEV of the energy-momentum tensor obeys 
the continuity equation {T^) c - U = 0, which for the geometry under consideration takes the form 

r^-(T 1 1 ) c + 2((T 1 1 ) c -(T 2 2 ) c )=0. (58) 

In addition, for a massless spinor field this part is traceless and the trace anomaly is contained 
in the point-like monopole part only. 

The fermionic condensate can be found from the formula for the VEV of the energy- 
momentum tensor by making use of the relation Tff = Mipip. The condensate is presented 
in the form of the sum of idealized point-like monopole and core-induced parts: 

(o|v^|o> = (o m |v^|o m > + (59) 

By taking into account formula (1531) for the components of the energy- momentum tensor, for 
the part coming from the non-trivial core structure one finds 



Wr^ J M xx ^ R ^ {ax) 



w)c = ^-2-irE z / dxx E 



.(60) 



This formula may also be derived directly from the mode-sum formula (0|^|0) = } 
with the eigenfunctions (|21|) . In deriving formulae (|53|) and (|60|) we have assumed that no bound 
states exist. In Appendix [A] we show that these formulae are valid also in the case when the 
bound states are present. 

For r > a the core-induced parts in the VEVs of the energy-momentum tensor and the 
fermionic condensate, given by Eqs. (|53p and (160f) . are finite and the renormalization is necessary 
for the point-like monopole parts only. Of course, we could expect this result as in the region 
r > a the local geometry is the same in both models and, hence, the divergences are the same 
as well. 

As it has been already mentioned, if there is no surface energy-momentum tensor on the 
bounding surface r = a, then the radial derivatives of the metric are continuous and, hence, 
in Eqs. (|32[) and (|57p one has u' a = 0, h' a = 1/a. In models with an additional infinitely 
thin spherical shell located at r = a, these quantities are related to the components of the 
corresponding surface energy- momentum tensor t^ u . Denoting by n^ 1 the normal to the shell 
normalized by the condition n^n^ = — 1 and assuming that it points into the bulk on both sides, 
from the Israel matching conditions one has 

{K^-Kh^} = 8itGt^. (61) 

In this formula the curly brackets denote summation over each side of the shell, hup = 9^v+n^n v 
is the induced metric on the shell, = h^h^V p n$ its extrinsic curvature and K = Kj^. For 
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i' a = 8nG U - ^ , K = - a + 4^Gr ° 



the region r ^ a one has = 5j l e v ^ and the non-zero components of the extrinsic curvature 
are given by the formulae 

K$ = - u '(r)e- v{r \ Kl = iff = -ti(r) e - v(r \ r = a-. (62) 

The corresponding expressions for the region r ^ a are obtained by taking in these formulae 
u(r) = v (r) = 0, h(r) = ln(ar) and changing the signs for the components of the extrinsic 
curvature tensor. Now from the matching conditions (I6ip we find 

(63) 

Note that the combination in the square brackets in Eqs. (|32p and (|57p is related to the surface 
energy-momentum tensor by the formula 

au' a /2 + ati a - 1/2 = 2irGar + 1/2 . (64) 

where r is the trace of the surface energy-momentum tensor. 

4 Vacuum polarization in the flower-pot model 

In this section, as an application of the general results given above we consider a simple example 
of the core model assuming that the spacetime inside it is flat. The corresponding model for the 
cosmic string core was considered in Refs. |23[ [24"1 [25] and for the global monopole core in Ref. 
|20j . Following these papers we will refer to this model as flower-pot model. Taking in the region 
inside the core u(r) = v(r) = 0, from the zero curvature condition one finds e h ^ = r + const. 
The value of the constant is found from the continuity condition for the function h{r) at the 
boundary which gives const = (a — l)a. Hence, in the flower-pot model the interior line element 
has the form 

ds 2 = dt 2 - dr 2 - [r + (a - l)a} 2 (d9 2 + sin 2 6d<j) 2 ). (65) 

In terms of the radial coordinate r the origin is located at r = tq = (1 — a)a. From the matching 
conditions (j63[) we find the corresponding surface energy-momentum tensor with the non-zero 
components and the trace given by 

r ° = 2r| = 2r| = ,~ \ 2vrGar = - - 1. (66) 

The corresponding surface energy density is positive for the global monopole with a < 1. We 
will consider the VEVs in the exterior and interior regions separately. 

4.1 Exterior region 

In the region inside the core the radial eigenfunctions regular at the origin are the functions 

Ri,nAr,k) = c/-^-, R 2 ,nAr,k) = c / X ° +n ° m , (67) 
yr yr 

where 

X a = j + 1/2 - n CT /2, (68) 
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and r = r + (a — l)a is the standard Minkowskian radial coordinate, ^ f ^ aa. Note that for 
an interior Minkowskian observer the radius of the core is aa. The normalization coefficient Cn 
is found from the condition 



M + uj 2kJ x ^(kaa) 



' P ^ouj Jl(ka) + Yl(kay 
with the tilted notation for the cylindrical functions 



F Ua {x) = xF' (x) 



' J' x (ax) 1/1 

X ~r — 

J\ a {ax) 2 



1 



a 



(69) 



(70) 



Note that Jy a {x) = for a = 1. In this case for the barred notation in Eq. (|57p one has 
f(~ 1,0 \x) = f( 1,0 \x). Hence, in the flower-pot model the part in the energy-momentum tensor 
due the non-trivial structure of the core is given by the formula 



oo 



ir 2 a 2 r 



5 V 



i=i 



dx 



M 



M 2 



C {ll-i/2(aax),I l/a _ 1/2 (ax)} ( ) 
G),i i/2 [ x i K l/a-i/2\xr)\ , 



C {l,_ 1/2 (aax), K l/a _ 1/2 (ax)} V**' 1 / 
where we have introduced the notation 



C{f(ax),g(x)} = xf(ax)g'(x) 



1 (\ 



a 



1 I f(ax) + xf'(ax) 



9{ x ) 



In formula ft]}, G#° [x,f u {y)} = F^' Vs) [xj v (y)\ for \i = 1,2,3, and 



^Z/i-l/2 [ X > fl 



fM\ = ( ^2 



1 



fl/a-1/2 (y) - ft, (y) 



Note that in terms of notation (|72p one has 

^W^W = ^^(ax),^^)}. 



(71) 



(72) 



(73) 



(74) 



For a = 1 we have C |/;_i/2(crax), Ii/ a -i/2{ ax )} = and the VEVs vanish. In the similar way, 
for the fermionic condensate from Eq. (|60p we find 



M 



oo 



7r 2 a 2 r ^— ' 



= 

M Vx 2 - M 2 



C {li-i/2(aax),I l/a _ 1/2 (ax)} - 



C{I, 



1-1/2 



'aax),Ki 



,X-l/a-l/2\ ax 



)} 



K l/a-l/2( xr ) ~ K l/ a +l/2( xr ) 



(75) 



For a massless fermionic field the core- induced part in the condensate vanishes. From formulae 
()7ip and (|75p it can be seen that for fixed values of r and M, in the limit a — > the core-induced 
parts vanish as a 2//Q+1 . Note that by using the recurrence relations for the modified Bessel 
functions, the functions C {f(oty),g(y)} in these formulae can also be written in the form 

c {ii-x/2(<*y), fi/a-x/2(y)} = ^fyii-i/2(ay)fi/ a +i/2(y) - 2/^+1/2(^)^/0-1/2(2/), (76) 

with f = I, K and A/ = 1, Xk = — 1- In particular, from formula (|76p with / = X it follows 
that C {-^-1/2(03/)) ^Q/a-i/2(2/)} < 0- As we W1U show in Appendix lAl this means that there 
are no bound states in the flower-pot model. 
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The core-induced part in the VEVs of the energy-momentum tensor and the fermionic con- 
densate given by formulae (|7ip and (|75p . are finite for r > a and diverge on the core boundary. 
Surface divergences are well-known in quantum field theory with boundaries and are investigated 
for various boundary geometries and fields. To find the corresponding asymptotic behavior for 
the points near the boundary, we note that in this region the main contribution comes from 
large values of I. By using the uniform asymptotic expansions for the modified Bessel functions 
for large values of the order (see, for instance, [26J), to the leading order we find 

(T °) c « -2(1%) c *—(Tl) c * 2 a ~/ (77) 
r — a 12Uir z aa{r — a) s 

,- v (l-a)M . , 

12-K z aa(r — a) 

As the parts corresponding to the geometry of the point-like global monopole are finite at r = a, 
from here we conclude that near the core the VEVs are dominated by the core-induced parts. 

At large distances from the core, r ^> a, in the case of a massless fermionic field we introduce 
a new integration variable y = xr and expand the integrand over a/r. The main contribution 
comes from the lowest mode I = 1 and to the leading order one has 

{T » )c ~ "SvrW (2 + a)I*(l/a + l/2) J Q VV L^^V-i/aWj ■ ( 79 ) 

Note that for a massless field the integrand does not depend on r. The integrals in Eq. ([79]) 
may be evaluated by using the formula for the integrals involving the product of the MacDonald 
functions (see Ref. [27]). For a massive field, assuming Mr ^> 1, the main contribution into 
the integral over x comes from the lower limit of the integration. Replacing the functions 
Ki/ a ±u2(xr) by the corresponding asymptotic formulae for large values of the argument, to the 
leading order we obtain the following estimates 

o ! ^J_/ T 2x „ exp(-2Mr) " ;2 C{/ M/2 HM),/ iM/2 (aM)} 

° \ l/c ~M/ 2/c ~ 4v / ^a 3 r 4 ^ C{l l _ l/2 {aaM) 1 K l/a _ l/ ^aM)}' 1 } 

As we see, in this limit the core-induced energy density and the radial stress are suppressed by 
the factor Mr with respect to the azimuthal stress. 

For a C 1 the solid angle for the exterior geometry is small and the corresponding scalar 
curvature is large. In this limit we replace the modified Bessel function containing in the index 
I jot by the corresponding uniform asymptotic expansions for large values of the order and the 
functions containing in the argument aax by the expansions for small values of the argument. 
In this way it can be seen that the main contribution into the VEVs comes from the mode I = 1 
and the VEVs are suppressed by the factor exp[— (2/ct) ln(r/a)]. 

In figure [1] we have plotted the dependence of the core-induced energy density (full curves) 
and radial stress (dashed curves) for a massless fermionic field as functions on the scaled radial 
coordinate r/a for a = 0.5 and a = 2. The azimuthal stresses are found from the zero trace 
condition. In figure [2] the same quantities evaluated for r/a = 1.5 are presented as functions on 
the parameter a. 

4.2 Interior region 

Now let us consider the vacuum polarization effects inside the core for the flower-pot model. The 
corresponding eigenfunctions have the form given by Eq. (|24p where the functions Ri, ntT ( r ; k) 
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1.2 1.4 1.6 1.8 2 



Figure 1: The core-induced energy density (full curves), a 4 (TQ) c , and radial stress (dashed 
curves), a 4 (T^) c , for a massless fermionic field in the exterior region of the flower-pot model as 
functions on r/a. The numbers near the curves correspond to the values of the parameter a. 
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Figure 2: The core-induced energy density (full curves), o 4 (Tq) c , and radial stress (dashed 
curves), a 4 (T 1 1 ) c , for a massless fermionic field evaluated at r/a = 1.5 as functions on the 
parameter a. 
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and R2,n rT (f, k) are denned by formulae (j67H . Substituting the eigenfunctions into the mode-sum 
formula, for the corresponding energy-momentum tensor one finds (no summation over ji) 



j=l/2 (7=0,1 U 



dx- 



J 2 (ax) J*(x) + YMx) 



(81) 



To find the renormalized VEV of the energy-momentum tensor we need to evaluate the difference 
between the VEV given by Eq. (|81|) and the corresponding VEV for the Minkowski bulk: 



<0i™ - (0mK\0 m ). 



(82) 



The appropriate form for the Minkowskian part is obtained from Eq. (|39p taking a = 1, replacing 
the integrand by [x, J\ a (xf /a)] and r — > r. As a result for the subtracted VEV one finds 



2tt 3 



aa 3 r 



i=i/2 
J^ 2 (ax)/a 



7T 



(83) 



The integral in this formula is slowly convergent and the integrand is highly oscillatory. 

In order to transform the expression for the subtracted VEV of the energy-momentum tensor 
into more convenient form, we note that the following identity takes place 



J x 2 (ax) I a 
J 2 /x)+Y 2 (x) 



7T" 7T" x - 

~4~ ~~ «~ 2-^ 



C{H{ s \ax),Hi:>(x)} 



8 s % C{J K {ax),Hl s }{x)} ' 



(84) 



Substituting Eq. (|84p into formula (|83p . we rotate the integration contour in the complex plane 
x by the angle ir/2 for s = 1 and by the angle —tt/2 for s = 2. Under the condition f < aa the 
contributions from the semicircles with the radius tending to infinity vanish. The integrals over 
the segments (0, iMa) and (0, —iMa) of the imaginary axis cancel out and after introducing 
the modified Bessel functions we can see that different parities give the same contribution. 
Consequently, the renormalized VEV for the energy-momentum tensor can be presented in the 
form (no summation over jj) 



00 



ir 2 af ' 
1=1 



dx 



x 



M 



M 2 



CjKt.yziaax),^^ (ax)} {fl) 

x ~Fut 7 7 — vT ( - T Z-i/2 L X ' 1 l-l/2{ xr )\ 

C{Ii-i/ 2 (aax),K l/a _ 1/2 (ax)} 1 l l l 1 ' 



(85) 



where the functions [x, f u (y)] are the same as in Eq. (|7Tj) with Xj = 1 in Eqs. ([55]) and (f56l) . 
Note that similar to Eq. (|76p. the function in the numerator of the integrand is also presented 
in the form 

C {Kl-l/2(<xy),Ki/ a -i/ 2 (y)} = -yK l _ 1/2 (ay)K l / a+l / 2 (y) + yK l+1/2 (ay)K l/a _ 1/2 (y). (86) 

For a = 1 one has C{Ki_i/ 2 (aax), Ki/ a _i/ 2 (ax)} = and, as we could expect, the VEV of 
the energy-momentum tensor vanishes. It can be explicitly checked that the components of the 
tensor given by formula (|85p satisfy the continuity equation (I58p and this tensor is traceless for a 
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massless field. In the way similar to that for the exterior region, for the renormalized fermionic 
condensate inside the core we find 



(V'V')ren = -5— T V" I / dx 



1=1 

C{K l _ 1/2 (aax),K l/a _ 1/2 (ax)} 



sjx 2 - M 2 

If_ 1/2 (xf)-/f +1/2 (xf)l . (87) 



C{Ii-i/ 2 (uax),K l/a __ 1/2 (ax)} 

The VEVs given by formulae (|85|) and (|87p are finite for f < aa and diverge on the core 
boundary. For the points near the boundary the main contribution comes from large I and we 
replace the modified Bessel functions by the corresponding uniform expansions for large values 
of the order. In this way it can be seen that the leading terms in the asymptotic expansion with 
respect to the distance from the boundary are given by the formulae 

/t>0\ o/t-^v 2aa irpiy a — 1 



aa — r 1207r 2 a 2 a(cra — f) 3 ' 

<**>- « vJj^-rV (89) 

Comparing these results with the corresponding formulae for the exterior region, we see that 
near the core boundary the energy density and azimuthal stresses have the same signs inside 
and outside the core, whereas the radial stresses have opposite signs. 

Near the core center the contribution of the summand with a given I behaves like f 2 ^ - ^ and 
the main contribution comes from the lowest mode I = 1 with the leading term (no summation 
over fi) 

(T u, „ r, * 3 C^(x) C{K 1/2 (ax),K 1/a _ 1/2 (x)} 

{ ^ /ren ~ 3ttW J Ma Vx 2 - M 2 a 2 C{I 1/2 (ax), K 1/a _ 1/2 (x)} ' 1 > 

where we have introduced notations 

C<°>(x) = 3(M 2 a 2 /x 2 - 1), C«(x) = C^ 2 \x) = 1. (91) 



In the case of a massless field, formula (|90p also gives the behavior of the vacuum energy- 
momentum tensor in the limit when the core radius is large and f is fixed, a/f > 1. In the 
same limit, for a massive field, assuming aM 3> 1, we replace in Eq. (|85p the modified Bessel 
functions containing in the argument ax by the corresponding asymptotic expressions for large 
values of the argument. By taking into account that the main contribution into the integral 
comes from the lower limit of the integration, to the leading order we have (no summation over 

-i —2aaM 00 

f 1 1 00 
( ^ )rcn W i^W^^^ 6 " 2 ^^^^^^-^^)], 03) 

/=1 

with f 7^ 0. In this case the energy density is suppressed with respect to the vacuum stresses 
by an additional factor (aM) -1 . 

For small values of the parameter a assuming that the core radius aa for an internal 
Minkowskian observer is fixed, we replace the functions Ki/ a _i/ 2 (ax) in Eq. (|85p by the cor- 
responding uniform asymptotic expansion for large values of the order. The leading term is 
obtained by making use of the replacements 



Cif^iaax), K 1/a _ 1/2 (ax)} -» yPn/,_ 1/2 (y)]' + ^l 2 + y 2 + 1/2, (94) 
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in the integrand of Eq. (|85[) with y = aax and f = I,K. As a result, in this limit the VEVs 
behave like 1/a. Due to the factor a 2 in the volume element the corresponding global quantities, 
such as total energy vanish as a. 

In figure [3] we have plotted the dependence of the renormalized vacuum energy density (full 
curves) and radial stress (dashed curves) for a massless fermionic field as functions on f/aa for 
a = 0.5 and a = 2. In figure H] the same quantities evaluated for f/aa = 0.5 are presented as 
functions on the parameter a. 
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Figure 3: The renormalized vacuum energy density (full curves), (cra) 4 (TQ ) su b> an d radial stress 
(dashed curves), (aa) 4 (T 1 1 ) su b, for a massless fermionic field inside the core of the flower-pot 
model as functions on f/aa. The numbers near the curves correspond to the values of the 
parameter a. 
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Figure 4: The energy density (full curves), (aa) 4 (TQ) su b, and radial stress (dashed curves), 
su b, for a massless fermionic field evaluated at f j aa — 0.5 as functions on the parameter 

a. 
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5 Conclusion 



In the present paper we have considered the polarization of the fermionic vacuum by the gravita- 
tional field of the global monopole wit a non-trivial core structure. The previous investigations 
in this direction are concerned with the idealized point-like model, where the curvature has 
singularity at the origin. In a realistic point of view, the global monopole has a characteris- 
tic core radius determined by the symmetry braking energy scale. For a general spherically 
symmetric static model of the core with finite thickness, we have evaluated the VEVs of the 
energy-momentum tensor and the fermionic condensate for a massive spinor field. These quan- 
tities are among the most important characteristics of the vacuum properties, which carry an 
information about the core structure. In the region outside the core we have presented the VEVs 
as a sum of two contributions. The first one corresponds to the geometry of a point-like global 
monopole and the second one is induced by the non-trivial structure of the monopole core. In 
the general spherically symmetric static model for the core, we have derived closed analytic 
expressions for the core-induced parts given by formula (|53[) for the energy-momentum tensor 
and by formula (|60p for the fermionic condensate, where the properties of the core are codified 
by the coefficient in the square brackets in the notation (I57j) . In Appendix [A] we show that the 
formulae for the core-induced parts in the VEVs of the energy-momentum tensor and fermionic 
condensate are valid also in the case when bound states are present. For points away from the 
core boundary these parts are finite and the renormalization is reduced to that for the point-like 
monopole geometry. Of course, we could expect this result as in the model under consideration 
the exterior geometry is the same as that for the point-like global monopole. For the points on 
the boundary the VEVs contain surface divergences well-known in quantum field theory with 
boundaries. 

As an example of the application of the general results, in section 0] we have considered a 
simple model with a flat spacetime inside the core. The corresponding model for the cosmic 
string core is known in literature as flower-pot model and here we use the same terminology for 
the global monopole. To have matching between the exterior and interior metrics, in this model 
we need the surface energy-momentum tensor located on the boundary of the core and having 
components given by Eq. (|66p . In the model with solid angle deficit (a < 1) the corresponding 
surface energy density is positive. The core- induced parts of the exterior VEVs in the flower- 
pot model are obtained from the general results taking as the interior radial functions in the 
eigenmodes the functions (|67[) . These parts are given by formula (|7ip for the energy- momentum 
tensor and by formula (|75h for the fermionic condensate. We have investigated the core-induced 
parts in various asymptotic regions of the parameters. In the limit when the core radius tends 
to zero, a — > 0, for fixed values r and M, these parts behave like a 2 ^ a+1 . For points near 
the core boundary the leading terms in the corresponding asymptotic expansions are given by 
formulae (|77p . (|78p . In this region the total VEVs are dominated by the core-induced parts. At 
large distances from the core these parts tend to zero as (a/r) 2//Qf+5 for a massless field and are 
exponentially suppressed by the factor exp(— Mr) for a massive field. We have also investigated 
the limit of strong gravitational fields corresponding to small values of the parameter a. In this 
limit the main contribution into the VEVs comes from the lowest mode I = 1 and the VEVs are 
suppressed by the factor exp[— (2/a) ln(r/a)]. 

For the flower-pot model we have also investigated the VEVs of the energy-momentum tensor 
and the fermionic condensate inside the core. Though the corresponding spacetime geometry is 
Monkowskian, the non-trivial topology of the exterior region induces vacuum polarization effects 
in this region as well. The renormalization is achieved by the subtraction from the mode-sums 
the corresponding quantities for the Minkowski spacetime. By making use of identity (|84|) . 
after an appropriate deformation of the integration contour, we have presented the renormalized 
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VEV of the energy-momentum tensor in the form (|85p and the fermionic condensate in the form 
()87p . These quantities are finite for strictly interior points and diverge on the boundary of the 
core with the leading divergences given by formulae (|88|) . (|89|) . In particular, near the core 
boundary the energy density and azimuthal stresses have the same signs inside and outside the 
core, whereas the radial stresses have opposite signs. Near the core center the main contribution 
comes from the mode I = 1 and the VEVs tend to a finite limiting value with isotropic vacuum 
stresses. Although the exact behavior for the fermionic field is unknown for a realistic model 
of the global monopole spacetime, the flower-pot model considered here presents some expected 
results as, for example, finite vacuum polarization effects at the monopole's center. For large 
values of the core radius the renormalized VEV of the energy-momentum tensor inside the core 
vanishes as a -4 for a massless field and as e~ 2aaM for a massive one. In the limit a <C 1, 
assuming that the core radius era for an internal Minkowskian observer is fixed, the vacuum 
densities in the interior region behave as 1/a. 

Note that in this paper we have considered quantum vacuum effects in a prescribed back- 
ground, i.e. the gravitational back-reaction of quantum effects is not taken into account. This 
back-reaction could have important effects on the dynamical evolution of the bulk model. We 
do not consider this extension of the theory, but note that the results presented here constitute 
the starting point for such investigations. 
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A Bound states 



In this appendix we consider the changes in the procedure described in the main text when 
bound states are present. For these states the quantity k is purely imaginary, k = i^y, and the 
corresponding exterior eigenfunctions have the form 

,,(e X ) _ n . ( . K ^Mn jlam (e, <p) 



with u> 2 = M 2 — 7 2 . To have a stable ground state we will assume that 7 < M. From the 
continuity of the eigenfunctions at r = a one has 

\/a-Ri,n CT (a, 27) = C h K Ua (-fa), ■ s /aR 2)riry (a,i^) = -in a C h K u<7+n<7 (>ya). (96) 

Excluding from these relations the normalization coefficient Ct> we see that for possible bound 
states 7 is a solutions of the equation 

KiyH ia )=0, (97) 

with the barred notation from Eq. (|57p . 

The coefficient Cb in Eq. (|95p is found from the normalization condition (|33|) . To derive a 
formula for the normalization integral, we rewrite equations (|16p . (1170 in terms of the functions 
Fp{ r ) = e u / 2+h fpir) and Gp{r) = e u / 2+h gp(r) and differentiate both equations with respect to 
u. Further we multiply the first equation by Gpfr), the second one by —Fp(r) and add them. 
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Combining the resulting equation with Eqs. (|16|) . (|17|) . it can be seen that the following relation 
takes place 



{Fj + G%) = e 



U+2/i 



9F, 



/3 



dG 



duj 



duj 



Integrating this relation we obtain the formula 



/ dre 

Jro 



v+2h\ f 2 



a u+2h 



M r )—g^ M r ) 



duj 



(98) 



(99) 



where ro is the value of the radial coordinate corresponding to the center of the core. By using 
this formula, in the case of bound states for the normalization integral in Eq. (|33p one finds 



i+2h 



( ,9Mr) .dgp(r) 

M r )—^ — M r ) 



[ dVT#V/J= I™ dr e v+2h [fl{r) + gl{r)} = e" 

J Jr L "~ J r =a+ 

(100) 

The expression on the right can be further simplified by using the continuity of the eigenfunctions 
on the boundary of the core. In this way for the normalization coefficient one finds the formula 

, 2 M + u 



n 2 



a- 



'<K v ^a)(d/du)K^(ja) 



(101) 



As a result, for the contribution of the bound state with k = 27 to the VEV of the energy- 
momentum tensor we have the formula (no summation over //) 



cc>, 



5 V 

7< 



■E'E 



7 



where v a = l/a — n a /2 and we have introduced the notations 



B^[ 7 ,K u ^r)], 



(102) 



b(%,k„m\ 



i 



M 2 
7*" 



M 



■v/M 2 - 7 2 



+ 



M 



VM 2 - 7 2 



+ 1 if 



v<T+n, 



i *i(y) 



BM[y,K Vir (y)] 



2ln a 



(103) 



ay 



KuAy) K u,+nAv), M = 2,3. 



In deriving Eq. (|102j) we have used the relation K Va {^a) = l//^' cr ' ) (7a) , which directly follows 
from the Wronskian relation for the modified Bessel functions in combination with Eq. (|97p . In 
the case when several bound states are present the sum of their separate contributions should 
be taken. 

The VEV of the energy-momentum tensor is the sum of the part coming from the modes with 
real k given by Eq. (|45D and of the part coming from the bound states given by Eq. (|102j) . In 
order to transform the first part we again rotate the integration contour in Eq. f|45[) by the angle 
7r/2 for s = 1 and by the angle —tt/2 for s = 2. But now we should take into account that the 
integrand has poles at x = ±^70 which are zeroes of the functions H^J (e Vs7T ' l / 2 x) in accordance 
with Eq. (1971) . Rotating the integration contour we will assume that the pole e nsm ^ 2 ^a, s = 1, 2, 

(s) 

on the imaginary axis is avoided by the semicircle C p in the right half plane with small radius 
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p and with the center at this pole. The integration over these semicircles will give an additional 
contribution 

—^hr y (2j+i) y y / dx^P-f^ M s) w«)L (104) 

j=l/2 ff=0,ls=l,2' yO f -"^ W 

By evaluating the integrals in this formula it can be seen that this term cancels the contribution 
(|102p coming from the corresponding bound state. Hence, we conclude that the formulae given 
above for the core-induced parts in the VEVs are valid in the case of the presence of bound 
states as well. 

In the flower-pot model the equation (|97p for the bound states takes the form 

C{l,_i /2 (aa7),«i/a-va(a7)} = 0. (105) 

In Section H] we have shown that the function on the left of this equation is always negative and, 
hence, in the flower-pot model no bound states exist. 
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